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We reveal that slow-light enhanced optical forces between side-coupled photonic-crystal nanowire waveguides can be flexibly 
controlled by introducing a relative longitudinal shift. We predict that close to the photonic band-edge, where the group velocity 
is reduced, the transverse force can be tuned from repulse to attractive, and the force is suppressed for a particular shift value. 
Additionally the shift leads to symmetry breaking that can facilitate longitudinal forces acting on the waveguides, in contrast to 
unshifted structures where such forces vanish. © 2011 Optical Society of America 

OCIS codes: (230.7370) Waveguides; (200.4880) Optomechanics; (350.4238) Nanophotonics and photonic crystals 



Development of reconfigurable photonic circuits utilizing 
gradient optical forces between micro- and nano-scale optical 
waveguides opens new possibilities for optical signal shaping 
and routing based on all-optical tuning of the structure geom- 
etry [T- "31. Optomechanical interactions can be enhanced in 
photonic-crystal waveguides due to increased optical energy 
concentration in the slow-light regime |[4| or through light 
trapping in cavities |5|. Whereas the pump intensity deter- 
mines the strength of gradient force, the sign of the gradient 
force between two identical coupled waveguides or cavities 
is dependant on the modal symmetry, with even modes typi- 
cally associated with attractive and odd modes with repulsive 
forces IS). This has been demonstrated in side-coupled wave- 
guides, where the force can be tuned from attractive to repul- 
sive by varying the phase of the incident light to excite modes 
of different symmetries |[Tl[2l. 

In this work, we suggest and demonstrate through 3D nu- 
merical simulations that slow-light enhanced optically in- 
duced forces between side-coupled photonic-crystal wire 
waveguides can be tuned by introducing a relative longitu- 
dinal shift. This happens because the shift breaks the struc- 
ture symmetry, which leads to pronounced hybridization of 
odd and even modes at the resonant frequencies close to the 
photonic band-edge where the group velocity is reduced. This 
enables the tuning of the transverse force from repulse to at- 
tractive, and force suppression for a particular shift value. We 
also reveal that the symmetry breaking can facilitate longitu- 
dinal forces acting on the waveguides, in contrast to unshifted 
structures where such forces vanish. 

We analyze suspended wire waveguides containing a peri- 
odic array of rectangular holes as illustrated in Fig. [Ha), and 
investigate the effect of the longitudinal shift zVx. To provide 
a link with experimental platform, we consider silicon nitride 
material (refractive index « ~ 2) and choose the dimensions 
according to Ref. ||5|, but without a cavity since we aim to 
utilize slow-light enhancement of interactions which can be 
realized in defect-free photonic-crystal waveguides l?). The 
structure sketch is presented in Fig. [T] The period of the rect- 
angular air hole array is a. The width and height of both wave- 



Fig. 1. (a) 3D sketch of the longitudinally shifted side- 
coupled photonic-crystal waveguides, (b) top view of the 
structure. 



guides are w = 7fl /6 and h = 2a/3, respectively, and the space 
between them \sd = a/6. The width and length of the rectan- 
gular air holes are Wx = 0.445a and w,, = 2a/3 respectively. 

We use a plane-wave expansion method 13 to calculate 
the supermode dispersion and the full 3D vectorial electro- 
magnetic fields. In Figs.|2a-c), we plot the dispersion curves 
for the first two transverse-magnetic (TM) bands close to the 
photonic band-edge, where the second band is the one we 
consider in the remainder of this paper. The corresponding 
group velocities for the second band are shown in Figs. IHd- 
f). The plots are presented in dimensionless units: normal- 
ized frequency (Oa/2Kc, normalized group velocity Vg/c, and 
normalized wavenumber ka/2K, where o is optical angu- 
lar frequency, Vg is the group velocity, k is the wavevec- 
tor component along the waveguides, and c is the speed of 
light in vacuum. The results are shown for three values of 
the longitudinal shift, as indicated in the figure caption. The 
group velocity vanishes at the photonic band-edge for all shift 
values, and this is a generally known feature. Additionally, it 
is also widely accepted that light-matter interactions can be 
increased in the slow-light regime, including an enhancement 
of the optical gradient forces j?)- However, as we demonstrate 
below, the introduction of a longitudinal shift has a nontriv- 
ial effect on slow-light enhanced optical forces, which can 
be attributed to breaking of symmetry of band-edge modes 
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Fig. 2. (Color online) TM mode properties in coupled wave- 
guides with different longitudinal shifts: (a,d,g) Ax = 0, 
(b,e,h) Ax = 0.15a, and (c,f,i) zix = 0.5fl. (a,b,c) Disper- 
sion relations shown as normalized frequency vs. wavenum- 
ber for the first- (dashed) and second-lowest (solid line) TM 
bands near the band-edge; grey shading marks the band- 
gap. (d,e,f) Group velocities for the second-lowest band. 
(g,h,i) Mode profiles at the band-edge: intensity (top 
row) and phase (bottom row), white contours show an out- 
line of the waveguides. 



from odd (anti-symmetric) profiles in the absence of a shift 
[Figs.[2g)], to complex phase profiles for shifted waveguides 
[Figs.|2l;h,i)]. 

We use the profiles of the optical modes to calculate the 
time-averaged gradient forces acting on the waveguides. We 
take into account the structure periodicity and associated 
Bloch symmetry properties of electromagnetic modes to sim- 
plify the general expressions based on the Maxwell stress ten- 
sor formulation r„jg (O. After excluding terms with vanishing 
contribution, we obtain the following expressions for the dis- 
tributed forces acting on the right waveguide in the transverse 
(y-axis) and longitudinal (x-axis) directions. 



fy - - 



dzJyy\y=(), /x = - 



dzr„\x='d- (1) 



In these expressions, the coordinate axis origin y = lies 
in the middle between the nanowire waveguides, and at this 
cross-section the tensor components are defined as 



r,,,|,=o = |£,|'-|£.p-|£.p- 



xy | v=0 — ^x^y + ^v^v + ^x^v + ^xi^y 



\H,\^-m^-\H,f, (2) 



(3) 



The force in the vertical direction vanishes, f- = 0, if the 
nanowires are free-standing or suspended high above a sub- 
strate. The gradient forces for the left waveguide have oppo- 
site values to those for the right waveguide given in Eq. ([T]i- 
We also analyze the net forces over one period, calculated as 



ajl 



fx,ydx. 
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Fig. 3. The force in the transverse (+y) direction on the right 
waveguide, (a) The force per unit length per unit energy den- 
sity vs. the longitudinal shift at the band-edge, (b-g) Forces 
for different longitudinal shifts: (b,c) Ax = 0, (d,e) Ax = 
0.15fl and (f,g) Ax = 0.5a. (b,d,f) The distributed force along 
the waveguide per unit energy density at the band-edge. 
(c,e,g) Total force per unit length and unit power vs. the 
wavenumber near the photonic band-edge. 



For comparison between different structures, it is conve- 
nient to consider a normalized dimensionless force [[6|, F = 
Fc/P, where P is the power flow through the waveguides. 
In the absence of losses the optical energy stored within one 
period of the structure can be expressed as U ~ Pa/\vg\. As 
the group velocity is reduced close to the photonic band- 
edge, the value of U increases, and the force value also grows 
as F ~ h'i'l order to compare forces in the slow-light 

regime for different structure parameters, we also consider a 
scaled force value F\vg\/c = F\vg\/P. 

In the transverse direction, the optically induced force orig- 
inates from the strong variation of the electromagnetic field 
due to the evanescent coupling between two waveguides. 
Most importantly, we show in Fig. 13 a) that the slow-light en- 
hanced force at the photonic band edge changes from repul- 
sive (positive value) to attractive (negative value) as the longi- 
tudinal shift is tuned from zero {Ax = 0) to half a period (Ax = 
0.5a). In the absence of a longitudinal shift (Ax = 0), the 
band-edge mode has odd symmetry about y-axis [Fig. I2g)] 
which provides a repulsive force at all waveguide locations 
[Fig- Kb)] in agreement with earlier studies for non-periodic 
waveguides |[6], and in periodic waveguides the force is addi- 
tionally enhanced close to photonic band-edge, see Fig.|3jc). 
When a longitudinal shift is introduced, the left-right sym- 
metry of the coupled waveguide structure is broken, and this 
leads to cross-coupling between odd and even supermodes. 
Consequently, the electromagnetic fields become a mixture 
of even and odd modes [Fig.|2h,i)] and the distributed force 
along the x direction changes its offset and peak-to-peak 
value, c.f. Figs. [3jb,d,f). As the longitudinal shift increases, 
the relative fraction of even supermodes gets larger, which 
explains the transition to transversal attraction between wave- 
guides at larger shifts. In particular, the strongest slow-light 
enhanced attraction occurs for half a period shift (Ax = 0.5a), 
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Fig. 4. The force in the longitudinal (+x) direction on the 
right waveguide, (a) The force per unit length per unit en- 
ergy density vs. the longitudinal shift at the band-edge, 
(b-g) Forces for different longitudinal shifts: (b,c) zVx = 0, 
(d,e) hx = 0.15a and (f,g) ^x = 0.5fl. (b,d,f) The distributed 
force along the waveguide per unit energy density at the band- 
edge. (c,e,g) Total force per unit length and unit power vs. the 
wavenumber near the photonic band-edge. 



when the distributed force is attractive at all waveguide loca- 
tions, see Figs.|3jf,g). Quite interestingly, for an intermediate 
shift value of Isx ~ O.lSa, the transverse force is suppressed 
even in the slow-light regime, see Fig. [3je), which happens 
because the distributed force oscillates between repulsive and 
attractive along the waveguide [Fig.[3ld)] providing zero av- 
erage force along a structure period. 

We estimate the characteristic values of the transverse op- 
tical forces, considering 100 mW optical power coupled into 
the waveguides and 1550 nm band-edge wavelength. Then, 
the transverse force per unit length is 4.45 nN//im for Ai = 0, 
a ~ 617nm, ka/2n = 0.499, and vjc = 0.0081 [Fig. He)]; 
0.02 nN/^m for Ai = 0.15a, a ~ 612nm, ;ta/27r = 0.499, and 
Vg/c = 0.0057 [Fig.He)]; and -39.91 nN/^um for /ir = 0.5a, 
a ~ 596nm, ka/ln = 0.4705, and vjc = 0.0018 [Fig.[3lg)]. 

In the longitudinal direction, the optical force also results 
from the strong electromagnetic field variation along the y 
direction according to Eq. The scaled force per unit 

length per unit energy density is shown in Fig.Ua). We note 
that it has a maximum value at A.^ ~ 0.3a shift, while zero and 
half-period shifts correspond to zero force. In order to explain 
this force dependence, we use the general symmetry proper- 
ties of electromagnetic fields in periodic photonic structures. 

For zero or half-period shifts, the refractive in- 
dex profile has a mirror symmetry with respect to 
a reflection transformation x — !> — x, and then it fol- 
lows m that Bloch wave profiles should satisfy 
relations Ey,{x,y,z;k) ~ Ey,-{—x,y,z\—k)e.yjp{i^\), 
Hy.z{^-,y-,z;k) = -Hy^,{-x,y,z;-k)sxp{i(pi), E^{x,y,z;k) = 
-Ejc{-x,y,z;-k)exp{i(pi), and Hy{x,y,z;k) = 
Hx{—x,y,z',—k)exp{i(pi). Additionally, for lossless 
structures with real refractive index profile, there is 
a symmetry related to time-reversal invariance ||9l: 
{E,H}*y,{x,y,z;k) = {E,H}j,^y,,{x,y,z;-k)exp{i(f)2). In 



these expressions the phase parameters (pi 2 can take or 
7Z values depending on the mode symmetry. Using these 
symmetry relations, the tensor component responsible 
for the longitudinal force according to Eqs. ([Hi,® can 
be transformed as: T^y = [—Ex{—x,y,z',k)Ey{x,y,z',k) + 
Ejt{x,y,z;k)Ey{~x,y,z;k) + Hx{-x,y,z;k)Hy{x,y,z;k) - 
H:,{x,y,z;k)Hy{-x,y,z;k)]exp{i(p) where 9 = (pi + (p2. 
It is then clear that Txy{x,y,z',k) = —Txy{~x,y,z',k), and 
accordingly the distributed longitudinal force is also an 
antisymmetric function of .ic, fx{x;k) ~ ~fx{—x;k). Indeed, 
the antisymmetric property of the distributed force is clearly 
visible in Figs. |4tb,f). Therefore, for structures with zero 
or half-period shifts the integral force over one period will 
always vanish, = 0, see the force dependence near the 
band-edge in Figs.|4jc,g). 

For intermediate shift values, < Ax < 0.5a, the structure's 
mirror symmetry is broken, resulting in a non-vanishing lon- 
gitudinal force, see Figs. |4td,e). For 100 mW optical power 
and 1550 nm band-edge wavelength, we estimate the force 
per unit length is \6.16yiN / pLm for Ax: = 0.15a, a ~ 612nm, 
ka/2K = 0.499 and vjc = 0.0057. 

In conclusion, we reveal that slow-light enhanced opti- 
cal forces between side-coupled photonic-crystal waveguides 
can be controlled by a longitudinal shift. The shift modifies 
the overall structure symmetry, enabling force control beyond 
the possibilities of symmetric structures with unshifted wave- 
guides. In the transverse direction, the longitudinal shift can 
be used to tune the optical gradient force from repulsive to 
attractive, and the force can be also suppressed. The shift 
also gives rise to force in the longitudinal direction along the 
waveguide, whereas this force component is always absent 
in symmetric structures. These results demonstrate that lon- 
gitudinal waveguide shifting presents new opportunities for 
controlling gradient forces which may contribute to the de- 
velopment of optomechanical circuits. 
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